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We calculate the early universe evolution of perturbations in the dark matter energy density in 
the context of simple dark sector models containing a GeV scale light mediator. We consider the 
case that the mediator is long-lived, with lifetime up to a second, and before decaying it temporarily 
dominates the energy density of the universe. We show that for primordial perturbations that enter 
the horizon around this period, the interplay between linear growth during matter domination and 
collisional damping can generically lead to a sharp peak in the spectrum of dark matter density per¬ 
turbation. As a result, the population of the smallest DM halos gets enhanced. Possible implications 
of this scenario are discussed. 


I. INTRODUCTION 

An important question for the field of particle cosmol¬ 
ogy is to reveal the nature of dark matter (DM). There 
has been compelling evidence for the existence of DM, 
most of which are from gravitational effects. There are 
also various ongoing experiments aiming to discover the 
DM through possible interactions other than gravity. On 
the theory side, there are well-motivated frameworks like 
supersymmetry that can accommodate DM candidates 
as weakly interacting particles. An alternative approach 
to the theory of DM is using minimality as guiding prin¬ 
ciple. The focus of this work is a class of simple models 
in which the DM candidate talks to the Standard Model 
(SM) sector through the exchange of a light mediator. 
This setup has been proposed and explored for several 
different motivations [IHS]. With a light mediator it is 
possible to suppress the direct detection constraints on 
DM, but still have prospects for other indirect searches. 
The simple interaction between the light mediator and 
DM can offer rich dynamics in the dark sector [THini- 

In this note, we investigate the light mediator phe¬ 
nomenology and cosmology in the context of two simple 
models, where the mediator is either a vector or scalar 
boson. Its interaction with the SM sector is via photon 
kinetic mixing in the vector mediator case, and the Higgs 
portal in the scalar case. We are interested in the mass 
range of the mediator at GeV scale, and lifetime as long 
as a second. This is often a less studied region of parame¬ 
ter space because the portal interaction is so weak that no 
current experiment is sensitive to it. The aim of this work 
is to fill this gap by studying the consequences of the me¬ 
diator’s longevity and possible testable implications. We 
note that the mediator is long-lived because of the tiny 
portal coupling between the dark and SM sectors—this 
also implies that the two sectors never reach equilibrium, 
but evolve with their own temperatures. There could 
be an interesting impact on the cosmology, if the light 
mediator becomes non-relativistic and temporarily dom¬ 
inates the energy density of the universe before it decays 
and releases entropy. We calculate the evolution of den¬ 


sity perturbations of DM during the mediator dominance 
and the imprint of this period on the primordial power 
spectrum. Because the mediator has become matter-like 
when it begins to dominate the universe, the DM density 
perturbation could grow linearly. A competing effect is 
the collisional damping from the DM-mediator scatter¬ 
ing. We find that for perturbations entering the horizon 
during the mediator domination era, the linear growth 
effect wins over damping. We show a generic predic¬ 
tion of such long-lived mediator model is a sharp peak in 
the transfer function right before the exponential cutoff. 
This corresponds to a characteristic length scale of order 
10“'^ pc, that is, size of the minimal DM halos. The peak 
implies these mini halos are to be more populated than 
in the usual cold DM case. 

In the next section, we study the simplest dark sector 
models with a light vector/scalar mediator. We summa¬ 
rize the current experimental constraints on the light me¬ 
diator mass and its portal coupling to SM. The allowed 
region where the mediator can be long-lived and domi¬ 
nate the universe is highlighted. In section III, we do the 
calculation on the evolution of density perturbations in 
the early universe, paying attention to the modes enter¬ 
ing horizon around the light mediator dominated era. We 
obtain the relative transfer function by comparing our re¬ 
sults with those of the usual cold DM. In section IV, we 
estimate the characteristic mass and size of the minimal 
DM halos and discuss possible testable predictions. 
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FIG. 1. History of dark sector and perturbation growth. 
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II. SIMPLE DARK SECTOR MODELS 
WITH A LIGHT MEDIATOR 


To set the stage, we sketch in Fig. the major events 
happening to the dark sector and their time scales rele¬ 
vant for this study. 


Vector mediator case 

The first dark sector model we examine contains a mas¬ 
sive vector mediator V^, 

C = CsM + - igxV^)x - m^XX 

. ( 1 ) 

The DM held x is a Dirac fermion and couples to 
the DM vector current. Both are SM gauge singlets. A 
kinetic mixing (the k term) between V and the usual 
photon is introduced. 

In the early universe, the DM x and the mediator 
are in thermal equilibrium through the gauge interaction 
g^- We focus on the case of light mediator, ^ my- 
In this case, the DM freeze out is controlled by the an¬ 
nihilation XX VV. To give the correct relic density 
to symmetric DM, the thermal annihilation cross section 
should satisfy 


anni — 


3 X 10“^®cm^/s 
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( 2 ) 


where S is the DM number density dilution factor due to 
any entropy production after the freeze out. For asym¬ 
metric DM, the relic density is given the primordial asym¬ 
metry in X and x number, and the annihilation cross 
section has to be larger than the thermal value. 

After the DM freeze out, the light mediator takes 
over all the entropy in the dark sector. We assume the 
dark sector has its own temperature T', and define its ra¬ 
tio to the visible photon temperature, ^ = T'/Ty, at very 
high temperature. After DM annihilation, the comoving 
number is 


nv{Tj) 


ume g*s{Ted) 3 

67r2 g*s{rnt) 


( 3 ) 


where g*s{T) counts the relativistic SM degree of free¬ 
dom contributing to entropy density. For this assump¬ 
tion to be valid, the interactions between and the SM 
particles must be decoupled throughout the history of 
universe. Because the (kinetic) mixing k is the only por¬ 
tal that communicate the two sectors, it is sufficient to 
require that the inverse decay rate // —?► is less than 
the Hubble parameter at temperature of universe equal 
to my- We find for < 10 GeV it requires k < 10“^. 


^ Although the other scattering processes such a,s qq Vg may 
have larger rates at T > ray because F ~ T, the relevant quan¬ 
tity here is the ratio to the Hubble rate, T/H. Because H ~ T^, 
these constraints are always less important than inverse decay. 
A key ingredient for this to be valid is that the photon has a 
universal coupling constant to all fermions. This is not true for 
the Higgs portal discussed in the next model. 


As shown below, this is well satisfied if the light mediator 
is long-lived enough to dominate the universe. 

There are further phenomenological constraints on n 
and my. In order to be consistent with standard big- 
bang nucleosynthesis, a commonly used requirement is 
to make decay before 1 second or so. This gives a 
lower bound on k. El We use the numerical result on the 
decay rate in m including the hadronic channels. On 
the other hand, k also controls the direct detection cross 
section of DM x, 


CTSI — 




mz 


( 4 ) 


where = g^/(47r). For weak scale DM, the current 
LUX [14] experiment places the strongest upper limit on 
K. In order to satisfy both constraints, the mediator V/^ 
should be heavier than ~ 50 MeV. El This is shown in the 
panel (A) of Fig.[^ We also notice that with a weak scale 
DM, the direct detection constraint on K — my parameter 
space is stronger than many of the limits on models with 
dark photon only. On the other hand, if the DM is lighter 
than around 4 GeV, the direct detection limits from LUX 
and GDMSlite [TB] (which covers the lower mass region 
down to = 4 GeV) are both evaded. However, as 
shown in panel (G) of Fig. the other constraints from, 
e.g., Supernova 1987a, are also very strong [T7HT^ (for 
an exception see [5D])— not much parameter space could 
be opened in the lower my region. 

We want to highlight a less constrained region where 
10”^^ ^ K < 10“® and my is at GeV scale. In this 
case, the V particle is long-lived enough and becomes 
non-relativistic before it decays; the total energy density 
of the universe could experience a temporary period of 
U-matter domination. To find the minimum lifetime of 
Vf^ for this to happen, one can compare its energy density 
with that of the visible radiation species. 


py{Tj) = myny{Tj)e 


Pr{Tj) 


30 
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( 5 ) 


Here g*{T) counts the relativistic SM degree of freedom 
contributing to energy density. After T' < my, the me¬ 
diator becomes non-relativistic, and py scales as T~^, 
while for radiation pr scales as T~'^ up to the changes in 


^ If there is an additional very light scalar or fermion in the dark 
sector which the mediator can decay into, the constraint on k 
could be evaded. This could lead to extra contribution to effec¬ 
tive number of neutrinos, Aeff |12l I13| . In this work, we focus 
on the minimal model and do not consider such extensions. 

^ The momentum exchange for DM scattering on xenon nucleus is 
roughly 50 MeV. Thus for mediator lighter than this value, the 
cross section will get saturated by the momentum exchange, and 
the bound could be slightly weakened m- We find this only 
slightly affects the lower bound on the mediator mass. 










3 


(A) m^=lTeV 



(B) m^=lTeV 



mv (GeV) 
(C) m^=2GeV 


(GeV) 



0.1 

mv (GeV) 



(GeV) 


FIG. 2. Experimental constraints on dark sector models with light mediator, shown in the parameter space of photon/Higgs 
portal versus mediator mass. We want to highlight the red region where the mediator could dominate the universe before it 
decays. 


g^fs- In order for pv to be ever larger than pr, we find 
the lower bound on the lifetime, 

Tv > 0.004^"® ( iGev ) ‘ 

Note 0.004 second corresponds to around a temperature 
of the universe at 15 MeV. 


This is a p-wave annihilation and T^d ~ 'm^l2Q is the 
chemical decoupling temperature. Again we assume the 
dark sector has its own temperature T'. After DM freeze 
out the comoving (j) number is given by 


_ 9C(3)g" 9*siT,d) ^3 

27r2 g^simt) ' 
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Scalar mediator case 

The Lagrangian for a dark sector with scalar mediator is 

^ = ^SM + xi^X - m^XX + + VxXX^P 

- vy2)cj) + , (7) 


where the DM x is a Dirac fermion, and the mediator cj) 
is a real scalar, and both are SM singlets. 

In early universe x nnd (j) are in thermal equilibrium 
via their Yukawa coupling y^. For ^ m^, the DM 
relic density is controlled by the process XX 4>4’- For 
symmetric DM, correct thermal relic density requires the 
annihilation cross section satisfy 


_ Q-Ko^Tcd 

2m^ 


3 X 10-26cm3/s 


( 8 ) 


In this model, there are two ways that cj) could commu¬ 
nicate to SM. First, if the coefficient of cross quartic term 
X(j,h is large, 4> will experience a similar freeze out via s- 
channel Higgs exchange. Second, the term in the 
Lagrangian induces a mixing between cj) and the Higgs 
boson and is responsible for (j) to scatter and decay. Be¬ 
cause of the Higgs portal, the (^-top quark scattering is 
most important in bringing (j) into equilibrium with SM 
than inverse decay. As said, we are primarily interested 
in a similar parameter space as the vector mediator case, 
where X^h, P 4 >h are small enough so that they never bring 
(j) into equilibrium with the SM sector. This requires, 
^<ph ^ 10“^ and < 10“^, respectively. 

The BBN and direct detection could further constrain 
and m^j,. For the 4> decay rate we use a result includ¬ 
ing the hadronic channels mull]. The direct detection 
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cross section of DM x is 


CTSI — 




mjjn 


4 ’ 
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( 10 ) 


where = y^/(47r). It was noted in [6] that if DM is 
heavier than ~ 5GeV, (j) must weigh more than twice of 
the muon mass, > 210 MeV. This can be seen from 
panel (B) of Fig. 

Because there are fewer constraints on a light scalar 
mediator than a vector one, for DM mass at GeV scale, 
we find it possible to lift the direct detection constraint, 
and make (f) as light as MeV. There is still an upper bound 
on (()-Higgs mixing from the low energy K'^ —>■ 

—>• 7r+t'+£“, —>• processes [23], and the 

fix-target GHARM experiment [24], as shown in panel 
(D) of Fig. In together with the decay before BBN re¬ 
quirement, there is a narrow window still allowed for the 
parameter 10“^ ;< ;< 10“^. However, in this 

window, (j) will always be thermalized with SM particles 
at high temperature, which fixes ^ = 1. 

Similar to the vector mediator case, we are interested 
in the possibility that (j) temporarily dominates the total 
energy density of the universe before it decays. Using ^ 
we derive the minimal lifetime of (j) for this to happen, 

~(l&v) <“> 

This happens if 10“^° < < 10“®, and is high¬ 

lighted in the red region of Fig. (B). 

To close this section, we have summarized the exper¬ 
imental constraints on two simple models of dark sector 
with light mediator, V or (j), and showed the regions in the 
parameter space that allow the mediator to temporarily 
dominate the universe as matter. The next few sections 
are devoted to the study of the implication of this possi¬ 
ble era on the DM primordial spectrum. 


III. LIGHT MEDIATOR DOMINANCE 
AND PERTURBATIONS 

To calculate the perturbation to energy densities in 
early universe, we define the metric in the conformal 
Newton gauge 

ds^ = —(1 -I- + (1 + 2^)a^Sijdx''dx^ . (12) 

The perturbations satisfy dt = — $ in the absence of 
quadrupole moment. 

There are three species in the univers^ DM y, light 
mediator </> and SM radiation r. The radiation includes 
all relativistic SM particles which tightly couple to each 


^ In this section we will proceed the discussion a with scalar light 
mediator (j). The results for vector mediator case are very similar. 


m0=lOGeV, T^=isec, f=l 



FIG. 3. Evolution of the energy densities of DM y, light 
mediator 0 and SM radiation r. In the shaded region, the 
universe is (j) dominated. 


other. For each species, the stress-energy tensor and the 
perturbations are 


= {1 + w)p{l + 5)ui^u'^ + wp[l + , (13) 

where w = p/p, and = (1 —tk, v). The next-to-leading 
order quantities are 5 and v. 

Using the continuity of the decay of (j) into usual 
radiation with DM as a spectator can be described as [25] 


Px + ^Hpx = 0 , 

p^ + 3{1 + w^)Hp^ =-T^p^ , (14) 

Pr+dHp^ = r^p^ . 


where an over dot means d/dt. The Hubble expansion 
rate is 


^^2 SttG 

H — —— [p^ + P4, + Pr) ■ 


(15) 


We include the process where the light mediator (j) 
changes from radiation-like {w^ = 1/3) to matter-like 
(w^ —)■ 0) as the temperature drops. As assumed 
throughout this work, after DM freeze out the total num¬ 
ber of y no longer change, and so does that of (j) until it 
begins to decay to radiation species of the SM sector. 

We present the results using a set of sample input 
parameters, = 10 GeV = Isec, ^ = 1, and 
= 1 TeV. We also choose the initial temperature to 
be Tini = m^j26 ^ right after freeze out, and de¬ 
fine the Hubble radius at this time to be oq = 1. As a 
good approximation, = T^i at the time Oq. The 

initial conditions for (14) are 


2^2 


Pxio^i) = s X (5.4GeV) x T}b—g:^s{Tini)T^ 


-3 

ini 7 




.4 

ini 5 


Pr{cii) — 2Q5*(Uiinz)2/, 


i4 

ini 


where rp, = n^js = 0.9 x 10 is the baryon-to-entropy 
ratio today, and the dilution factor S due to entropy 
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ni^ = 10GeV, r^ = lsec, f=l m^ = 10GeV, T^ = lsec, f=l 



FIG. 4. Evolution of primordial perturbations of different wavelengths that enter the horizon before, during and after the (j) 
domination regime. Different colors correspond to wave numbers k/kf.q = 10®'® (cyan), 10®'® (red), 10®'® (magenta), 10® (or¬ 
ange), 10®^'® (yellow), 10®" (green) and 10® (blue), respectively. The dashed part of curves means sign change in the perturbation 
due to acoustic oscillation. 


release from </> decay is fixed a posteriori from solving 
Eqs. (fni). 


/ Px(ai) \ / Prjaf) 

VPx(«/)/ \Pr{a^) 


3/4 


(16) 


where af corresponds to an epoch well after (j) decay. 

Fig.j^shows the evolution of the p’s. During the period 
10^ ^ o, <3 X 10®, the energy density of (j) beats that of 
radiation and dominates the total energy in the universe. 
This happens until the </> particles decay away, and dump 
a significant amount entropy into radiation. With the 
above set of parameters, S ~ 60, and matter-radiation 
equality happens at aeq = 2.8 x 10 ®®. 

The equations for the evolution of perturbations in the 
fluid densities and in the metric are 


+ ^ + 3$ = 0 , 


(17) 


+ H9y H-$ — Cj. -dy +Tq^( 


-Ox) 


X) 1 




<^0 + (1 + ^4‘)~ , 

0^ + il- 3W^)H0^ + —$ - 

a l + w^ a l + w^ 

(5, + ^ + 4$ = {6^ - (5, - $) , 

6a Pr 

k'^ , 


Or - + —$ = 

4a a Pr 


$ 


H + 


e \ 


3Ha? 


$ = 


47rG 

w 


04, - 0J 


i.Px^X P(f>^(fi T Pr^r) 1 


where 0 = aV • v for each species. We derived these 
equations following the notations in [26U2^ . Notice again 
the over dot means d/dt here. The rate t“® is related to 
the scattering process (px —t (t>X which keeps the DM in 


kinetic equilibrium with p, and it is given by 



where the first bracket is the px px cross sectioij^ 
at T ^ TOy, and the second bracket represents the ra¬ 
tio of collision time to the thermal relaxation time |30) . 
The relative velocity = 1 when T > m^, and 
Vif, = y/3T/mrj) at low temperature. For = 0.1 and 
my = 1 TeV, this rate falls below Hubble at akd — 10^ 
(around T ^ lOMeV). With this input, the sound speed 
Cs of the the y fluid is, = 5p^/5p^ = 2T^/{3m^) be¬ 
fore kinetic decoupling (a < akd), and more suppressed 
for a > akd- The back reaction of the scattering on the 
00 equation has been neglected because p is much more 
populated than y. 

The primordial scalar perturbations does not evolve 
until entering the horizon. The initial conditions for 0 
can be obtained by simplifying these equations for super¬ 
horizon modes. At time oq, the universe is still radiation 
dominated, thus, 

3 3 

(5y(ao) = 1 ^ 0 ( 00 ) = -5r{ao) = -$o , 

Oxia-o) = 0^{ao) = 0 ^( 09 ) = > (19) 

where $0 is the primordial scalar perturbation in the 
metric. 

We solve the set of equations in 0 for modes that 
enter the horizon before, during and after the p domina¬ 
tion regime. The results on $ and dy are shown in Fig.|^ 


® For the vector mediator model, the Vx Vx cross section is 
that of Thomson scattering, Stto® /(3m® ). 
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FIG. 5. The transfer function normalized to that of collision¬ 
less cold DM. keq corresponds to the mode that enters the 
horizon at matter-radiation equality. 


The momentum exchange term ~ ^x) 

important because 6^ and 6^ evolve hand in hand with 
each other. 

To obtain the transfer function T{k), we evolve the DM 
density perturbation to the time of matter-radiation 
equality, and normalize the perturbations to those of 
a cold DM. The result is shown in Fig. The rela¬ 
tive transfer function T{k)/TcBMik) exhibits a peak at 
wavenumber ^ lO^keq right before the exponential cut¬ 
off. The enhancement factor in this case in the density 
perturbation can be as large as 100 — 1000. 

This peak-plus-cutoff feature is a unique prediction of 
dark sector models where DM couples to a long-lived light 
mediator. This result also does not depend on whether 
the DM is symmetric or asymmetric. 

IV. POSSIBLE TESTABLE IMPLICATIONS 


From the $ plot one could infer when a mode enters the 
horizon, while the plot shows the impact of (j)-X scat¬ 
tering and the era of (j) matter domination on different 
scales. 

During the period when > pr, the universe is tem¬ 
porarily matter dominated, and it is well known that in 
this period the density perturbations of (p itself and DM x 
both grow linearly. This feature is shown in Fig.[^ where 
(j) domination depicted by the shaded region. For modes 
entering the horizon in this period, <i> drops slowly allow¬ 
ing 6^ to grow, which then serves as the source for the 
linear growth of S^. After </> begins to decay and the uni¬ 
verse returns to radiation dominated, these perturbations 
grow logarithmically. Because all (p particles eventually 
decay away, only the DM density perturbation 6^ mat¬ 
ters for later structure formation at this scale, and it still 
remembers the (/)-domination history. The enhancement 
in 6^ implies it will reach the nonlinear growth regime 
much earlier than the usual cold DM case. 

At earlier time, the effect of (p-x scattering plays an 
important role. The scattering term couples the (p 
fluid tightly to the DM x fluid. Because in this regime (p 
still carries non-negligible kinetic energy and has a larger 
sound speed, it drives the acoustic oscillation and causes 
collisional damping in 6^ until the rate falls below H at 
akd- At very high wavenumber A:, the resulting damping 
factor in the 5^ spectrum can be estimated as [5^ . 


k^dampiPkp — exp 



( 20 ) 


where we And kcut corresponds to the mode entering 
the horizon at the moment of (p domination begins. For 
modes entering the horizon during p dominated era, they 
are less affected by the collisional damping effects, al¬ 
though the kinetic decoupling happens afterwards. This 
can be understood because both p and x have become 
effectively non-relativistic, and their perturbations sat¬ 
isfy the same equation with the same initial condition. 


From DM density perturbation spectrum one can cal¬ 
culate the DM structure formation at later time. The 
peak before the cutoff in the primordial spectrum cor¬ 
responds to a special scale of the smallest DM ha¬ 
los/clumps. The peak right before the cutoff indicates 
that these smallest objects are most likely formed. Here 
we estimate their properties. 

From Fig. we get the wavenumber at the peak in 
spectrum, kpeak ~ W^k^q. The mass of the smallest ob¬ 
jects today can be determined kpeak ^ evaluated using |31) . 

Mpeak = f ~ 10-^ Mq . (21) 

^ \ >^peak J 

where we used today’s DM relic density = 0.26, 
today’s critical density = 8.5 x 10“^°gcm“^, keq = 
0.073Mpc“^D5)M^^i k = 0.673, and Mq is the solar 
mass. 

The size of these halos is given by the length scale of the 
mode k~^^ at the time when the DM density perturbation 
Sx{o,ni) approaches to 1 [3^ . 


p _ I - 1 

Peafe \4TrPelm{anl) J 


( 3M 


peak 


1/3 


( 22 ) 


Clearly Rpeak is more sensitive to the linear growth dur¬ 
ing p domination. To get ani, we first assume scale in¬ 
variance in the primordial scalar perturbation and fix 
the initial value of $o — 8 x 10“® from the Harrison- 
Zel’Dovich-Peebles spectrum in the large scale structure 
observation [34]. From this, we find for the mode corre¬ 
sponding to the peak of Fig. S^ikpeak) will reach non¬ 
linear growth shortly after recombination, and Rpeak ~ 
10^"^ cm. The average DM density inside such mini halos 
is about 10® times higher than the local DM density in 
the usual case. 

Therefore, the first formed DM halos have a typical 
mass similar to that of the earth, and typical size similar 
to the earth-sun distance. After the formation, these ha¬ 
los may merge into larger objects or be destroyed by grav¬ 
itational tidal forces. A more precise calculation of the 










7 


small scale structure formation and properties of such ob¬ 
jects today calls for simulation and is beyond the scope of 
this work. Today, if the smallest halos still takes substan¬ 
tial fraction of total Milky Way galaxy mass, the average 
distance of two nearby such objects is ^ 10^® cm. Thus 
the solar system would run into such objects every 50- 
100 years and the direct detection rates will get enhanced 
periodically. Another way to observe these clumps is in¬ 
direct detection signals from, e.(?., DM annihilation from 
inside. For asymmetric DM case, if the self-interaction 
allows DM bound states to exist and form, as discussed in 
the context of Yukawa bound states [5] or the case of dark 
atom [^, the indirect detection signal is also available. 
It could be enhanced if the bound state formation makes 
these clumps dissipate enough and further reduce their 
sizes. There have also been suggestions that large pri¬ 
mordial density perturbations could lead to micro-lensing 
effect if DM forms more compact objects [36] . 


V. OUTLOOK 

To summarize, in the context of simple dark sector 
models containing a long-lived light mediator, we notice 
the mediator’s longevity can be closely connected to an 
era of temporary matter domination. We have calcu¬ 
lated the early universe evolution of dark matter den¬ 
sity perturbations, especially during the mediator dom¬ 
ination era when the perturbation could grow linearly. 
We showed this scenario can result in a sharp peak in 
the transfer function compared to that of collision-less 
cold dark matter. This leaves us with more structures 


on small scales to look for, and could provide additional 
handle in testing light mediator models with very weak 
interactions to the standard model. More dedicated sim¬ 
ulation of the small scale structure formation would help 
to obtain more concrete prediction in the properties of 
these mini dark matter halos. 

We note a temporary matter dominated universe be¬ 
fore BBN has been discussed in various warm DM mod¬ 
els, such as wino DM produced from late decaying mod¬ 
uli [37], and keV sterile neutrino DM diluted by long- 
lived right handed neutrinos [55HiD] . In the moduli decay 
case, the free streaming of boosted wino always washes 
out most structures on small scales [?7l [^ . The free 
streaming nature of light sterile neutrino will result in 
the same suppression. In contrast, in the light media¬ 
tor models discussed in the present work, the key for the 
primordial small scale structures to survive is that the 
long-lived mediator cannot decay to DM, and the DM 
itself remains cold. 
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